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Abstract 

We consider 2-dimensional QCD on a cylinder, where space is a circle of length L. We formulate 
the theory in terms of gauge-invariant gluon operators and multiple-winding meson (open string) 
operators. The meson bilocal operators satisfy a W^ current algebra. The gluon sector (closed 
strings) contains purely quantum mechanical degrees of freedom. The description of this sector in 
terms of non-relativistic fermions leads to a W^ algebra. The spectrum of excitations of the full 
theory is therefore organized according to two different algebras: a wedge subalgebra of W^ current 
algebra in the meson sector and a wedge subalgebra of W^ algebra in the glueball sector. In the 
large N limit the theory becomes semiclassical and an effective description for the gluon degrees of 
freedom can be obtained. We have solved the effective theory of the gluons in the small L limit. We 
get a glueball spectrum which coincides with the 'discrete states' of the (Euclidean) c = 1 string 
theory. We remark on the implications of these results for (a) QCD at finite temperature and (b) 
string theory. 



0. Introduction 



Soluble models of string theory are not abundant and this greatly limits our ability to tackle 
many of the important problems in string theory. In the last several years there has been consid- 
erable amount of work in the c = 1 string theory fij. However this theory, being a closed string 
theory in two dimensions, has only one propagating degree of freedom and hence does not provide 
an opportunity to understand the physics of the infinite tower of massive states which are in higher 
dimensions responsible for some of the most remarkable properties of string theory. Recently there 
has been a trend to explore lower dimensional gauge theories [g, y, |4|, ||. One motivation for this 
is to provide new models of string theory. In a previous paper Q we had considered the problem 
of QCD with fermions in a two dimensional spacetime which is a plane. We formulated the theory 
in terms of gauge-invariant open string operators (which physically represent a meson) satisfying a 
Woo algebra. We reinterpreted the infinite tower of mesons (a result originally due to 'tHooft @ ) as 
perturbations around a classical ground state that form a representation of the wedge subalgebra 
Woo+ <8> Woo + of the original Woo algebra, thus providing an example of an infinite dimensional 
algebra that organizes the massive levels of string theory. 

In this paper we discuss the more complicated example of QCD (with fermions) on a two- 
dimensional cylinder where space is a circle of length L and time is non-compact. Unlike the 
case of the two-dimensional plane, QCD2 on a cylinder has dynamical gauge fields, besides the 
Coulomb force between the quarks, giving rise to a much larger class of gauge-invariant operators. 

(n) 

Let us denote by Txy the path which starts from point x on the circle and reaches the point y 

after winding around the circle n times in an anticlockwise fashion. The non-abelian phase factor 

corresponding to this path is given by the operator U^ n '(x,y, t)\ a b = (P exp[ig J p („) Ai(x f ,t)dx']) a b- 

xy 
The present theory is described in terms of two basic sets of gauge-invariant operators constructed 

out of the above operator: (1) the pure gluon operators U(n,i) = jjtrU^ n >{x,x,t) and (2) the 

multiple-winding 'meson' operators M^ J/3 (n,t) = j^^ ab ^f a {x,t)U^ (x^y^)^^^) (here i, a 
and a respectively denote flavour, dirac and colur indices). Because of the presence of non-trivial 
winding numbers, the 'mesons' this time satisfy a Woo current algebra (W^) (Sec. 2) and obey 
constraints that depend on the U(n,t) (Sec. 3). These constraints define a coadjoint orbit of W^ 
and lead to a Kirillov construction of the classical action (Sec. 4). 

One of the interesting features in the quantum mechanics of our model is the nature of the 
fermi vacuum (Sec. 5). We find that there is no fixed fermi vacuum labelled by a fixed fermi level 
since the latter is neither a gauge-invariant notion nor a conserved quantity due to the anomaly 
in the axial current. The right vacuum instead is a linear combination of fermi vacua labelled by 
all possible values of the fermi level. The implication of this on the gluon dynamics (which, in the 
gauge A\ = constant and diagonal, is described in terms of N non-relativistic fermionic particles 
moving on a circle) is that effectively the fermions live on the real line under the influence of a 
harmonic oscillator potential and a mutual interaction (Sec. 6). Some of these phenomena have 
earlier been observed by Manton [H] in the context of the Schwinger model on the cylinder. 

The description of the gluon sector in terms of non-relativistic fermions leads to the interesting 
phenomenon that there are two Woo algebras in our theory; one is the current algebra W^ of the 
'mesons' described earlier and the second is the Woo algebra that is formed by fermion bilinears 
where the 'fermions' now mean those in the gauge sector. Indeed since these latter fermions live 
on a circle whose length goes as 1/L (as against the quarks which move in real space of radius 
~ L) there is a rather interesting duality between these two sets of fermions. In the small L limit 
the 'meson' fluctuations (excitation energies ~ !/-£) are effectively frozen (Sec. 8) and the Woo of 



the gluon sector becomes operative. We find the 'glueball spectrum' in this limit and show that it 
coincides with the 'discrete states' of the c = 1 string theory and is arranged by the W m algebra 
|]]. In the large L limit, on the other hand, it is the mesons which fluctuate and the gluons are 
frozen. In that limit the physics reduces to that of the planar model where the meson fluctuations 
are arranged by a W^ algebra || which is the zero mode part of the current algebra (Sec. 9). 
For intermediate values of L we indicate how meson fluctuations around the ground state form 
a representation of a wedge subalgebra of W^ (Sec. 7) although the interpretation of them as 
spacetime particles and how they interact with the glueballs remains an open question. 

1. The action and the hamiltonian 



We consider the gauge group U(N). As usual we denote the gauge fields, which are hermitian 
matrices, by A" where \x = 0, 1 is the Lorentz index and a,b = 1,2, ... ,N are colour indices. The 
fermions are denoted by ipf a where i = 1,2, ... ,nf is the flavour index and a = ±1 is the dirac 
index. We consider one space and one time dimension where the space dimension is a circle of 
length L. 

The theory is described by the following path-integral 

Z = J VA (x, t)VAi{x, t)Vip(x, t)V^(x, t) exp[iS(A , A x ,ip, ip*)) 
S = / T dt / L dx( - \trF^F^ + ^(iD^ - m^) (1) 

F 01 = E = d t A x - d x A + ig[A , A x ], D^ = <9 M + igA^ 
which has the following gauge invariance 

ip(x,t) -> ip n {x,t) = Q(x,t)ip(x,t) 

A^(x,t) -► A%(x,t) = id^{x,t)^{x,t) + 04^t 

where £l{x,t) is an arbitrary U(N) matrix. We first fix the gauge d x A\ = 0. Using the standard 
Fadeev- Popov procedure we can reduce VA\{x,t) to T>A(t) with appropriate Jacobian factors. We 
get 

Z = JVA{t)A 1 [A}VA (x,t)ViP(x,t)V^(x,t)exp[i(S + S coul + S F )\ 

AiL4] = det(id x /g + ad A) 

S = Jdt^[tT(d t A + l g[A^,A]) 2 } (3) 

Seoul = JdtJ^dxl[tv(D A A ) 2 - 2gtv(pA )], Pab {x,t) = ^ b {x,t)^ a {x,t) 

S F = Jdt J L dx[^id t ip + tp^ 5 (id x - gA)^ - mip^°ip] 

In the above we have denoted the constant mode of the field Ao(x,t) by Aq (t) and the rest of it 
by Ao(x,t). We still have residual gauge transformations which are purely time-dependent, Q(£), 
which we now use to further gauge fix Aq to zero. Going through the Fadeev-Popov procedure 
once again we get 

Z = /Pi(t)A 1 [i]PA)(x,t)2?V(x,i)PV t (^,t) exp[i(5 + S cout + S F )] 
S = J dt -§tr(c>M) 2 



alongwith the constraint 

Lp ,aa = Qa = (5) 

Here po ab are the generators of colour rotations on fermions and are zero modes of p ab (x) appearing 
in (H), thus 

Po>ab = (1/L) f dx^ b (xW{x) (6) 

Jo 

The Coulomb term S cou i and Sf remain the same as in (J3|) . To proceed further we use the following 
change of variables: 

A(t)=V(t)X(t)vHt) (7) 

where the colour matrix A(£) is diagonal. The measure T>A changes as 

VA(t) =PA(t)W(t)n t A 2 (A(t)) 

W(t) ^Y[ t Y[ a ^(t)dV(t)) ab (8) 

A(A(t)) = Ua<b(^ - h) 

We have used a constraint (V^dV) aa = to make the change of variable (^) well-defined. Under 
this the action Sq becomes 

So = f dt±[tr(d t \) 2 - Ea^(Aa - \ b )%bha} 

U = {v ] {t)id t v{t)) ab 

In order to compensate for the change of variable (J7|) in Sf and S cou i we make corresponding 
changes of integration variables ip — ► V'iJj, Aq — » VT^oV- The integration measures are invariant 
under this change. Sp and S cou i now become 

S F = J dt J Q L dx[^idt^ + tr(pol) + ^j 5 (id x - g\)^ - m^W] 
S co «z = / ^/ L dxiftrCDAA)) 2 - 2ptr(pIo)] 

Here we have defined p ab = p a b{x) — po.afe- The term tr(/?oO appears because of the change of 
fermion variable ip(x,t) — ► W(t)^(x, i). Note that Z aa does not appear in it because of eqn. (5). 
That is consistent with the change of variable (7). The functional integration over V(t) requires a 



careful treatment [11|. We get 



flit A 2 (A(i))W(t) exp N J J dt± J2 a ^b [(Aa - \b)Xbha + 2l abPoM ] 

= A(\(0))A(X(T))exp(iS ,coui) ( n ) 

So,coul = ~~ Ji o ^% Sa^fe P0,abP0,fea(A a — Aft) 

It is easy to do the integration over Ao(:c, t). This produces a Jacobian factor l/Ai[A(£)] which 
cancels the Jacobian in the measure in (H). Note that AiL4(i)] = Ai[A(i)]. In the action this inte- 
gration amounts to replacing in S cou i in (|10|) Ao(x,t) by —g(D\)~ 2 p. In terms of mode expansions 
tp(x,t) = J2n V'n(i) exp(2irinx / L) this means replacing in S cou i the classical value 

<^(2w5irf-W' " #0 (12) 



Combining the resulting S cou i with the off-diagonal zero mode contribution So tCOU i of (11) and 
collecting all the pieces we get 

Z = jVX(t)A(X(0))A(X(T))V^(x,t)V^(x,t)exp[i(S + S F + S coul )] 

S F = J T dt J L dx[ipUd t %b + ^h 5 {id x - g\)ifj - mift'fil)] 
Seoul = J2a,b Jo dt J L dx J L dyp ab {x)p ba (y)K ab (x - y) 

In the above the kernel K ab {x) is defined by 

K ab (x) = -g 2 L/(8<K 2 ) En=-oo exp(-i2imx/L)[n + (X a - X b )/X )} 2 , a^b 

(14) 

K aa {x) = -3 2 L/(87r 2 )En^oexp(-z'27rna;/L)/n 2 



where 



Ao = ^ (15) 

gL 



Note that the action in ( |i~3| ) is invariant under the "large" gauge transformations |1, 12 1 

fL ab (rh, x) = exp[-igxm a X ]5 ab 

(16) 

>% = ^a + m a X 

The measure is not invariant under these transformations because of the presence of the van der 
Monde factors A(A(0)) and A(A(T)). The transformations ( |J6[ ) were part of the residual gauge 
transformations in the gauge d x Ai = and were broken when we fixed the gauge A Q ' =0. 
This is because the mode expansion of Aq is not invariant under the large gauge transformations 
(16): Aq' — > AQ b n+rna _ mb . Now that Aq has been integrated out completely, we can restore this 
invariance by observing that since the classical action is invariant, the path integral automatically 
invariantizes the van der Monde factors (this can be established by an argument involving a change 
of integration variable corresponding to ([!(])). The result is that A(A(0)) and A(A(T)) are replaced 
by their periodic counterparts Ap(A(0)) and Ap(A(T)): 

A P (A) = nsm[vr(A a -A fe )/Ao] (17) 

a<b 

The fact that this is the right measure can be inferred in two steps. First, we know this to be the 
case for pure Yang- Mills theory in two dimensions [i"4| , |l5| ]. In the presence of dynamical fermions, 
since the fermion measure is separately invariant under large gauge transformations, and the new 
terms Sp and S cou i in the action are also invariant, the invariantization of the measure for the 
gauge fields must proceed in an identical fashion, giving rise to the same formula (fL7|). 

The Ap factors at t = and t = T in the measure ensure that the initial and final wavefunctions 
are completely antisymmetrized with respect to the A a 's (since to start with they must be symmetric 
with repsect to permutation of the A a 's on account of Weyl-symmetry) . This gives rise to the well- 
known fermionic nature of these eigenvalues. 

The kernels 



The kernels (14) can be explicitly evaluated. In the (x,y) region of integration pertinent to (13) 



the argument of the kernel K ab {x) lies in the range [— L, L\. In this range 

Kab(x) = -(g 2 /4)exp[igx(\ a - \ b )][L/ (2 sin 2 [it (\ a - \ b )/\ ]) - ixcot[ir(\ a - A 6 )/A ] - |x|], a^b 

K aa (x) = -(g 2 /A)[L/6-\x\+x 2 /L], x e [-L,L] 

(18) 



As we have observed earlier, the action in (13) is invariant under the large gauge transformations 
(|16|). We would like to rewrite the theory in a form that makes this symmetry manifest. Let us 
introduce the periodic (invariant under (Jig)) 5- function 



5p(X-X a )= J2 £(A-A a + nA ) (19) 

n=— oo 

and note that 



/ ° d\5p(X-X a ) = 1 (20) 

Jo 

Using this and the invariance of S cou i under ( |16| ) we may write K a ^ b (x) as 



K a ^ b {x) = f X ° dX [ X ° dX'K(X - A'; x)e^<^-^e i ^ x "-^5 P {X - X a )5 P (X' - X b ) (21) 

jo Jo 

where 

K{X;x) = -{g 2 /4)e i9xX [L/(2sm 2 [-irX/Xo]) - ixcot[irX/X ] - \x\] (22) 



Note that we cannot use the ctp-functions in the integrand in ( |2l| ) to set A a = A and A;, = A' in 
the exponential factors since these 5p-functions only equate the fractional parts (in units of Ao) of 
respectively A a and X b with A and A'. We use the uncancelled integer (in units of Ao) parts of A a 
and Aft in the following to write the contribution of (ET1) to S cou i as 



/ dt f dx f dy f ° dX f ° dX'K(X-X';x-y)Y,Pab(x,t)p ba (y,t)5p(X-X a )5p(X' - X b ) (23) 
Jo Jo Jo Jo Jo ~T, 



dX'K{X-X'-x-y) 

where p ab {x,t) = x^ b {x,t)x a (x,t) and x a ( x ,t) is defined as 

Xta(x,t) = exp[igx(X a - XM a (x,t) (24) 

Note that the x' s are invariant under the large gauge transformations (16). Now, the a ^ b sum 
in ( p3[ ) can be thought of as J2 a b ~ 2~2 a =b- The latter, when inserted in (f23|), involves the quantity 
K(X; x) in the limit A — > 0. In this limit K(X; x) behaves as 



K(X- )X )\ X ^ = - 9 - 



L/2 L . . x 



2 



+ 



(vrA/Ao) 2 '6 ' ' L 



+ o(X) (25) 



Because of the zero charge condition, eqn. (|J), the divergent term in ( p5| ) does not contribute to 
(|23]) . The rest of it exactly cancels the K aa (x) contribution to S cou i. Hence, we get 



Seoul = dt dx dy dX dX'K(X-X';x-y)]>p ab (x)p ba (y)5p(X-X a )5p(X'-X b ) 
Jo Jo Jo Jo Jo ~L 

a,b 

(26) 



The bilocal 'meson' operator 

Introducing now the manifestly gauge-invariant bilocal field 

^(A.*) = iE^H(A- W^^t) (27) 

a 

we may rewrite ( P6|) in the manifestly gauge-invariant form: 

Seoul = —^\ dt dx dy dX dX'K(X-X';x-y)tr df {M xy {X;t)M yx (X',t)) (28) 

4 JO JO JO JO JO 



where 

L/2 

sin 2 (-7rA/Ao 



K ( A ! x ) = 2/--X /x n ~ ^cot(vrA/A ) - \x\ (29) 



and the trace 'tr^/' is over dirac and flavour indices. We notice that in terms of the manifestly 
gauge invariant operator M(A, t), an internal dimension of length Ao has emerged. 

The physical meaning of the gauge-invariant bilocal fields introduced in (£7|) is transparent in 
the exponential representation of the periodic 5- function (Poisson's summation formula): 



oo 



5 P (X -X a ) = ^- V e i27m ( A - A ^ Ao (30) 

An 

u n=— oo 

Substituting this in (g7j) we see that 

-. oo 

M^(X,t) = Y J2 e-^ nX ^M^(n,t) (31) 

^ n=— oo 

where 

M^(n, i) = 1 X; e^^^C'C^ *)^(y, t) (32) 

a 

is the bilocal "meson" field with n windings of the gauge field around the spatial direction. We 
note that ( p2|) is indeed the manifestly gauge invariant operator 



M^in, t) = 1 X <,(*, «) ( e" ^ ^ (x '' i)dx ") ^t* (y> t) (33) 

a V /aft 

where Iiy is a path connecting a; and y with n windings around the circle^. In the gauge 

A 1 (x, t) = V\t) X{t)V(t), [X{t)] ab = S ab X a (t), (35) 



it simply evaluates to (|32|). We also note that the periodicity of the ^'s in the spatial direction 
implies the constraints 

M^in, t) = M^(n + l,t), M^{n, t) = M l ^(n - 1, t) (36) 



1 In higher dimensions T xy is not characterized by an integer but M xy (T,t) still form a closed algebra (see next 
section) : 

[M xy (T,t),M x , y ,(F',t)} = S xy ,M y , x {V'oT,t) - S x , y M xy ,(ToT' \t) (34) 

where ToT' is a curve defined by (ToT') xy i = F xy T' x , ,S yx i . 



In terms of the variable A these read as 

M^ J/3 (A, t) = e i27rX / Xo M^ j P(\, t), M$? (\, t) = e~ i27rX / x °M^ jP (X, t) (37) 

The hamiltonian 

We may now write down the Hamiltonian for this problem 

— = Hym + Hp, (38) 

where 

Hym = WlT,pI ( 39 ) 

a 

and 

H F = J X ° dXTr [ifdM{X) - m 7 °M(A)] - *£ J L dx J L dy J X ° d\ J X ° dX'K(X -X';x-y) 

ti df (M xy (X)M yx (X')) 



K(X;x) is given by fl29|) , d xy = d x 5p(x — y) and Tr stands for trace over dirac and flavour indices 
and integration over spatial coordinates. The periodic spatial <5p-function is defined similarly to 

oo 

6 P (x-y)= J^ S{x-y + nL) (41) 

n=— oo 

The above hamiltonian comes with the prescription that it acts on wavefunctions which already 
include the factors Ap(A) described in (O) and are therefore completely antisymmetrized in the 
Vs. 



2. Operator algebra of the bilocal field 



The bilocal fields (|27D (equivalently (p2|)) satisfy a closed algebra. This can be easily derived 
using the fermion anticommutation relation 

{iPt a (x,t)^(y,t)} = S ab S ij 6 a0 S P (x-y) (42) 

The gauge invariant operators (p2)((|33|)) satisfy the following current algerba 

[M™M(n, t),M l X/?\n\ t)] = (l/N)^' 5^' 5 P {y - x^M^/^'in + ri, t) / ^ 

, , (43) 

- (l/N)tiF 8°P 8p(x - y')M l / y ^(n + n', t) 

or equivalently the loop algebra 

[M X ^(X, t), M%J' ' (A', t)] = (l/N)^'S^'Sp(y - x')6 P (X - X')M x f^' (A, t) 
- {l/N)& ST? 8p(x - y')5 P (X - X')M^{X,t) 



(44) 



For n = n! = 0, we recognize (Q) as the infinite dimensional algebra W<4 = W^ ® U(nt) ® 
U(2). (E3) ((|43])) involves an infinite dimensional generalization of the algebra Woo, which may 



appropriately be called Woo current algebra. We shall denote this in the standard notation by W^ 
while W^q will stand for Woo <8> U{rif). We have seen in || how a central term arises in W<4 algebra 
of the bifocal meson fields (in the case in which L —* oo) when expanded around a given classical 
background. A similar phenomenon occurs in the present case and fluctuations of M(A) around 
any given classical background satisfy a loop algebra with a central term. This central term is 
essentially the ground state expectation value of the bifocal field M(X). This is discussed in detail 
in Sec. 7. 



3. Constraints 

We now show that the zero charge condition, eqn. (|5|), implies a quadratic constraint on M(A). 
Using fermion anticommunication relations we can prove the identity 

(M 2 (A))^ = (l/iV)5>p(A-A a )(M(A))^ (45) 

<8< 



+ (L/N 2 ) J2 e^-^ x 5 P (\ - \ a )5 P (\ - \ b )xUx)x%(y)P0,a» (46) 

a.b 

where po,ab, defined in eqn. (||), is the generator of colour rotations on fermions and x' s are the 
gauge-invariant fermions defined in (^J). Here we have used the convenient matrix notation for 
the square of M(A), £a, 7 / L dz (M(A))£'* 7 (M(\))%' i(3 = (M 2 (A))^. Now, because of the 5- 
functions in the second term in (|4q), the sum over a,b receives contribution only when A a = A^. 
However, since the physical states are antisymmetrized with respect to (A a , A&) for every pair (a, b), 
this can happen only for a = b. Hence, the double sum in ( ^6|) collapses to only the a = b terms on 
physical states. But this contribution also vanishes on physical states because of the zero charge 
condition, (g). We, thus, arrive at the constraint 

M "W = [jr E ^( A " V*))]M(A) (47) 

a 

Let us introduce the density of eigenvalues 

W(A,t) = ^E^(A-Aa(t)) (48) 

a 

By definition, it satisfies the constraint 

Ao 

dXU(X,t) = l (49) 

o 

To see the physical meaning of the eigenvalue density in the pure Yang-Mills sector, we write using 

-1 oo 

U(\,t) = — V e- i27TnX ^°U(n,t) (50) 

An — ' 

u n=— oo 

Now U{n,t) is the gauge-invariant loop variable with n windings in the gauge (|35"|): 

U( n , t) = —Y e «««A.(t)/Ao = 1 tre igf T (.n) Aidx (51) 

N ^-^ N 

a 

9 



r( ra ) is a curve that winds n times around the circle; — n corresponds to windings in the opposite 
direction. Using (^) in (|47| ) we finally arrive at the constraint 

M 2 (A) =U(\)M{\) (52) 



Also, from the definition (27) of M(X) and the zero charge condition, eqn. (|5|), it follows that 

Tr(M(A) - W(A)) = (53) 



(p3) and fl53| ) constitute an infinite set of constraints which the bilocal variable M(X) must satisfy. 
As we shall see later, these constraints specify a coadjoint orbit of the group W^ and enable us to 
construct a classical action in terms of the bosonized variable M(A). 

4. The large N limit — classical action and equations of motion 



In the rest of this paper, we will consider the present model in the limit iV — > oo. Since 
quantum fluctuations are suppressed in this limit, the theory greatly simplifies. The large iV limit 
is taken by holding g = gvN fixed as N — > oo. In this limit, then, Ao grows as yN . This is 
because Ao = 2-ir/gL = \/N 2n/gL = \JN Ao, where now we must hold Ao fixed as N — ► oo. This 
necessitates scaling of A, A a and p a to respectively X/^/N = A, X a /^J N = X a and p a /VN = p a . 
Note that in terms of the scaled variables the value of h is 1/A^. Henceforth, we shall drop the 
'bar' on Ao, A and A a . We shall, however, continue to use g = g\N to remind ourselves that in the 
following A^ — ► oo limit has been taken. 

We also define here another variable M(A) that we shall actually need in the following discussion: 



(m(a)Y qj/3 = e-^-y^ (M(x)) ia,jl3 

V / xy V / xy 

= (l/NU(X,t))Z a Sp(X-X a )xt a (x,t) X %(y,t) (54) 

= (1/NU(X, t)) E a 5 P (X - X a ) exp[ig(x - y)(X - A a M a (x, t)i/>%(y, t) 

where M(X) = M(X)/U(X). By construction, M(X) is periodic in x,y. It is, however, not periodic 
in A. 

In terms of M(A) the fermionic part of the hamiltonian, Hp, eqn. (|38|), is given in the limit 

A^ — > oo by 



pA„ 



H F = Jq° dX U{X) Tr j 5 (id - gX) M{X) - m 7 u M(A) 



■> 



~ L a„ r L a^ r x o a\ ij(\\ r x o 



~ 9 T So dx Jo d y So dX U W So dX> U(X') K(X -X';x- y) tr df (M xy (X)M yx (X')) 
where 



(55) 



K{X;x) = e iSxX K(X;x), (56) 



and K(X;x) has earlier been defined in (^). Also, in terms of M(X) the constraints (|52j) and ( p3|) 
read 

M 2 (A) = M(A), (57) 

Tr(l-M(A)) = 0. (58) 



10 



Classical Action and Equations of Motion: 

In order to find out the classical action of the system, we need to recognize the important fact that 

p a and M do not commute. Indeed they satisfy the commutation relation 

p a (t),M(X,t)] =-^d Xa M(\,t) (59) 



From the definition of M(X,t), eqn. (|54j), we see that the right hand side of (p9|) does not vanish. 
Had this commutator been zero, it would have been easy to write down independent symplectic 
forms for the (X a ,p a ) system (namely, simply J2a.dX a A dp a ) and for the M system (OTj) and the 
classical action would involve sum of these two symplectic forms. We shall see below that the 
(^a,Pa) system carries a representation of the Virasoro Group (Diff(S' 1 )) while the M system, as 
already noted, forms an orbit of Woo- The existence of the commutator (|59| ) owes to the fact 
that the Virasoro group acts on the M system as well, thus leading to an intertwining of the two 
group actions. The situation here is similar to the one encountered in rigid body dynamics where 
the space-fixed angular momenta and the linear momenta do not commute. To disentangle these 
two group actions we shall adopt below a procedure analogous to going to the body-fixed angular 
momenta. 

Virasoro Action : 

The Virasoro group arises as follows. Let us consider the space of the A a 's. The change from one 
set of A a 's to another set, A' a 's, can clearly be achieved by a reparametrization / : R — > R, f(X a ) = 
X' a Va. The easiest way to see this is to think of a two-dimensional plot in which the A a 's are 
marked on the x-axis and the X a : s are marked on the y-axis. The points (x, y) = (A a , A^) provide 
N number of points for the plot of y = f(x). In the limit N — » oo the map is completely specified. 
Let us now consider only those changes of A a 's in which the integer parts of A a do not change. In 
other words, if we introduce the useful notation 

x = int(x) + frac(x) (60) 

where int(x) is the largest integer less than or equal to x and frac(x) = x — int(x), we are looking 
for functions / which are such that int(/(A a )) = int(A a ) Va. It is clear that these functions are 
actually maps from S 1 — > S 1 and are of winding number one. These functions form a group (the 
Virasoro group) under standard composition and inversion. The reason for restricting to these 
functions is that in the discussion of equations of motion we will be concerned with infintesimal 
time-evolutions in which int(A a )'s do not change. 

It is easy to see that the density U{X) of eigenvalues which is a single-valued function on the 
circle A £ [0, Ao] changes under a map / _1 : S 1 —>■ S 1 , A a i— »• X' a = / -1 (A a ) as follows 

U f (X) = f'(X)U(f(X)) (61) 

where we have defined 

U f (X)^U(X\{f-\X a )}), U(X)=U(X\{X a }) (62) 

Thus a change of the eigenvalues A a — » / -1 (A a ) is equivalent to a reparametrization A — » /(A). 
Under these reparametrizations U{X) transforms as a spin-one current. It is convenient to introduce 
the quantity <p(X) defined by 

U(X) = d xV (X) (63) 
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(/9(A) transforms under the above Virasoro group as a scalar: 

<</ {f {\)) = <p(\) (64) 

We now make the important observation that the entire U(X)- configuration space is an orbit of the 
Virasoro group. The reason is that, as we have argued, any distribution A a of eigenvalues can be 
reached from a fiducial distribution A aj o by an appropriate S 1 — ► S 1 map g~ x , that is 

Aa = <?~ 1 (A„,o) (65) 

provided that the integer parts of these two distributions are the same. Let the densities cor- 
responding to the fiducial distribution and the arbitrary distribution be called Uq(X) and IA{X) 
respectively. By (|61|) we have 

U{\) = U {X)dg/dX (66) 

By choosing the fiducial distribution A aj o to be such that Uq is uniform (by the normalization 
condition it must then be 1/Ao), we can explicitly determine the Virasoro transformation g~ l that 
takes us from A aj o to A a . The result is 

dg/dX = XqU{\) (67) 

Using ( |63[ ) we find 

<?(A)/A = f X U(X) = ip(X)-ip(0) (68) 

Jo 

Thus the orbit of the Virasoro group can be parametrized by the values of the function <p(X). 
We can repeat the above exercise for M. We get 

M(/(A)|{A a }) = M(A|{r 1 (A Q )}) (69) 

where, as in fl62|) , we have displayed the dependence of M(X) on {A a } for clarity. Let us choose / to 
be precisely the transformation g^ 1 that takes us from the fiducial distribution A aj o to A a , namely 
(pq). In that case, we see that 

M( 5 - 1 (A)|{A a }) = M(A|{A ,a}) (70) 

In other words, if we define a new variable 

M(X) = M(g- 1 (X)) (71) 

then according to (f70|) .M(A) does not depend on the {A a }. More precisely, 

[ Pa ,M(X,t)}=0 (72) 

This is the disentanglement we were looking for. 

Because of the zero commutation relation ( f72| ) it is now easy to write down the "pq" term for 
the whole system as the sum of the separate "pg" terms for the {X a ,p a ) system and the A4 system. 
The new field M satisfies the same constraints as M 

M 2 (X)=M(X) (73) 

Tr(l-M(X)) = (74) 
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and the W<£, algebra 



An 

N 



5 P (\ - A') {s^S^'Spiy - x')M™/ P \\) - 5^5^' 5 P {x - y^M^^X)} 



(75) 



A kinetic term can be constructed for it using a coadjoint orbit construction, similar to that used 
in |1, Il0| , for the present case of W^ . In this way we arrive at the classical action 



S = N 



i / E dsdt Jq q dXM(X, t, s) [d t M(X, t, s), d s M(X, t, s)] 

+ I dt{ZaPadtX a - (E„ JlPI + Hp)} 



(76) 



where the region X of (s, t) integration is the lower half plane, £e(— oo, +oo), se(— oo, 0] and we have 
the boundary condition that M.(X,t,s = 0) = A4(A,t) and A4(A,£, s) — ► t-independent function of 
A as s — ► — oo. In the above equation Hp is given by (|55|). 

We can readily obtain the equations of motion from ([?(]) by the standard variational principle. 
The allowed (by the constraints (|73|) and (|74|) ) variations of ,M(A, t) are the W^ group rotations, 



M(X,t) -► e il4/ ( A '')^(A,t)e- il4/ ( A ' i ) 

= A4(A,i)+i[W(A,i),M(A,i)]+o(PF 2 ) 



(77) 



After obtaining the equations of motion it is convenient to go back to the original variable M(X,t) 
and reexpress these equations in terms of it. After these manipulations we find that (^6[) leads to 
the equations 



(dt-iEaPadx^If^^t) 



-7 5 (w 



K(X' -X;x-z) 



Id - gX) + m 7 °, M{X, tj\ ***** + f J A ° dX' U(X't) J L dz 

J xy 

M xz (X',t)M zy (X,t)) -K(X-X';y-z)[M xz (X,t)M zy (X',ty 

d t X a = -j-, d t pa = -d Xa H F , 



where the kernel K(X;x) is given by (pq) and (p9|) . For completeness we reproduce it below: 



K(X;x) 



igx\ 



L/2 



sin 



m 



,'ttA 

zxcot I — 
Ao 



(78) 



(79) 



(80) 



Note that the right hand side of the second equation in fl79f) is non-zero because of the implicit 
dependence of Hp on {A a } through U(X) and M(A). 



5. The Anomaly and the fermi vacuum 



In this section we would like to study the effective theory that describes the dynamics of the 
jauge fields obtained by 'integrating out' the fermions. For this purpose, in the large N limit that 
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we are considering, it is enough to restrict ourselves to the fermi vacuum. Actually as we shall see, 
there is really no single fermi vacuum in this theory in the conventional sense. One has to sum 
over fermi vacua labelled by all possible values of the fermi level. This is intimately connected with 
the presence of anomaly in the axial current in this theory as is explained in detail below. In the 
following we shall set the quark mass to zero and consider the case of a single flavour only. 

Let us recall that the fermionic part of the hamiltonian, given in (|55[), has a four-fermi interaction 
term. As we shall see later, this term can be neglected in the large N limit for the discussion of the 
effective gauge theory since the error thus introduced vanishes in this limit. The quadratic fermion 
part (linear in M) of (|55| ) can be rewritten in the massless case in terms of fermion modes as 

Hf,o = -£E T K7 5 <V4 a )(n + 5 A a ) (81) 

a n 

We note that the operators O^ = Tr{j WnVn) commute with Hpfl for all a and n. In order to 
construct the eigenfunctions of H we can therefore start by constructing simultaneous eigenfunc- 
tions of the 0® . Since we are interested in the ground state of Hp^ we need to consider only those 
eigenfunctions of the 0% which minimize Hf,q for any given background values of the A a . Now 
an eigenfunction of the O^ which minimizes Hpfi in a given background of A a 's is nothing but the 
filled fermi sea in that background. We therefore proceed to construct the latter. 

The fermi sea in a fixed Yang-Mills background: 

To construct the fermi sea we consider, as usual, the free dirac equation for quarks of each colour 

and dirac index 

[idt + J 5 (id x -g\aW a = (82) 

Choosing 7 = o\ and 7 1 = —i(J2, we have 7 5 = 7°7 1 = 03. In this representation the dirac 
equation becomes 

[id t + sga(a)(id x - g\ a )W a = (83) 

where a = ±1 represent the right-, left-moving fermions respectively. The above equation is solved 
by 

<(x, t) = £ expH^t + i^x\r a , n (84) 

n L 

where 

E a na = ^sgn(a)(n + ^) (85) 

describes the single-particle energy spectrum. For a = +1, E% n increases with n. Thus the fermi 
(dirac) sea is constructed by filling the levels from n = —00 upto some fermi level n° =1 = n^. For 
the left-moving fermions, similarly, the fermi sea is constructed by filling from n = +00 down to 
the fermi level n^ = _ 1 = n\. 

Fermi level transforms under large gauge trnasformations: 
We recall that under the large gauge transformations (16) 

r a (x, t) - {^)%{x, t) = expH^^)C(z, t) (86) 



which implies 

Ufl\a _ 

y a,n+m a 



(^L=C,n +m „ (87) 
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Thus the defining condition for the fermi level n R 

(^ + a i, n r + i, n )=0(-n + n a R ) 
leads to the following equation for the gauge-transformed fermionsfj 

((^+l,n^ U )+l,n) = (4V+m>+l,r,+mJ = 0(~n - m a + n R ) (89) 

From ( |88| ) and (^9|) we see that the fermi level (n R ) for the gauge-transformed fermions is given 
by 

(n R f = n R -m a (90) 

Similar statements are also valid for the fermi level for left-movers and one gets 

(nlf =n a L -m a (91) 

The left- and right-moving fermi levels are related 

We observe that the left- and right-moving fermi levels are not independent and they actually get 

related by the condition Q a = (see equation (||)). The value of Q a in the fermi sea is 

n R OO 

Qa = £ i + £ i ( 92 ) 

n=— oo n=n^ 

This is a divergent sum and it needs to be regularized. In the following we will regulate divergent 
sums over fermionic modes (of the colour index a) by the gauge-invariant regulator exp(— e\n + 
^a/Ao|) which cuts off contributions from large values of \n\ []. The result is 

Q a (e) = ~ + (n R + l-n a L ) (93) 

We will define the regularized charge Q„ e9 by subtracting the divergence, that is 

Q? a = n% + l-ni (94) 

This vanishes only if 

n a L =n a R + l\fa (95) 

The subtraction that we have used here and the resulting expression for the regularized charge is 
consistent with gauge invariance. There is a beautiful argument due to Manton [|| which shows 
that a violation of (p5l) leads to a non-gauge-invariant answer for the total momentum which is 
formally gauge-invariant on states of total charge zero. His argument is for the abelian (Schwinger) 
model but it can be generalized to our case quite easily. In the following we shall use the notation 

n a R = n a F , n a L = n a F + l (96) 



2 In ( |8q ) and (183) we have introduced theta functions with integer arguments, defined such that 8(n) = 1 for 
n > and otherwise. Note that with this definition we have the identity 9(n — 1) + 0(—n) = 1. 

3 This regulator gives results that are non-analytic in X a for finite e. However, the finite e-independent pieces of 
regularized quantities can be shown to be analytic in A a . 
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Expectation value of the quark bilinear in the fermi sea 

Let us define the notation |{n^}) for the fermi sea in a given background of A a . The state |{n^}) 

is defined by the following expectation values 

(r+i,n^i, n } = ^n-n a F -l) (97) 

and 

(fiii) = «h+4) (98) 



Using the above results it is easy to calculate the quantity M xy {\) (eqn. (|54|)) in the fermi vacuum 
\{n a F }). We get 



OO 



«■ — - ■ 27rn(x — y) 

M xy {X,t) = - Y, e J ^^M n (A,t) (99) 

n=— oo 

where 

M n (A) = j^ J2 E 6 ( X - & + mX °) K^^W" -m-l) + (l^fL)e(-n + m)\ (100) 
In the above equation we have defined a new gauge-invariant quantity 

ia = Aa + U F (101) 

. 27r(x — y) 

Note that because of non-periodicity of M xy (X,t) in A (M xy (X + Xo,t) = e l l M xy (X)), the 
modes M n {\) satisfy the following quasi-periodicity condition 

M n (A + A ) = M„+i(A) (102) 

The above condition is clearly satisfied by M n (X) calculated in (p_00|). 



It is easy to show that M n (X,t) obtained from ( 100 ) by replacing £ a by £ a (t) and U{\) by 



U(X,t) satisfies the constraints (|57[) and the (58), and the equation of motion (]78|). The quadratic 



constraint is satisfied quite trivially. The second (trace) constraint is satisfied because we have 
already ensured Q a = 0. To see that the equation of motion (J78|) is also satisfied, we note that 
M xy (\, t) in (|99l ) is a function of (x — y) only. This fact, together with the dirac structure of 
M n (A, t) displayed in ( |100| ), ensures that for massless quarks the right hand side of ( [78| ) vanishes. 
The fact that the non-linear term in M does not contribute at all to the equation of motion is 
the classical counterpart of the statement that the fermi sea constructed by ignoring the quartic 
fermion term is correct in the N — » oo limit. With the right hand side equal to zero, equation ([78]) 
reduces to 

(dt-yJ2Padx a )M n (X,t) = (103) 



L „ 



It can be easily checked that M n (\,t) obtained from (|100[) by replacing £ a by £ a (t) and U(X) by 
U(X,t) satisfies (ffOg). 



Effective hamiltonian for the gauge fields: 

The effective hamiltonian for the gauge fields is obtained from the full hamiltonian H by sub- 
stituting for M(X) its value (|100| ) in the fermi vacuum (note that the contribution of the excited 
states is subleading in 1/N). The pure Yang-Mills part of the hamiltonian of course remains the 
same. H F becomes (with quark mass set equal to zero) 



H F = gX f ° d\U(X)V(X) - — J™ / ° dA^(A) / ° d\'U(X')K{\, A') (104) 

Jo 2LA n Jo Jo 
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where 



and 



VW = ~ E (A + nA )tr d/ ( 7 5 M n (A)) (105) 



°^ tr d/ M„(A)M n /(A') 

*<*.*> = E [(A+n [ Aa) _ (A , +n , A) p doe) 

n,n'=— oo Lv ' v /J 

with M(A) given by ( |100| ). Using the exponential regularization introduced above we get 



2 1 
e ~ 12 



K(A) = - T - — + U re9 (A) + o(e) (107) 



and 



^ 2 



where 



K ^( A ' A ') = ■ 2 ,tm 77T + 2 ( ln e + C - 1) + iW A , A') + o(e In e) (108) 

esm A^i A- A J 

K eff (A) = ^^y E 5 P (X - CaXWAo + 1/2) 2 (109) 

K reg (X, A') = N 2k(x)k(\') 51 MA-£a)<Jp(A-&) W™a&) + V'(-^afe) + w ab {ip'(w ab ) - ip'{-w ab )}] 

(no) 

Here w ab = (£ a — £fc)/Ao and C is Euler's constant. Also, ip(x) = 4- lnT(x) and ij/(x) = -^ip(x), T(x) 
being the standard gamma- function. Once again, £ a = A a + n F , the gauge-invariant combination. 

Now notice that the above expressions depend on the fermi level n F . The dependence of 
quantities such as total energy on the fermi level is of course quite natural. However, since in 
our problem the n F 7 s are not gauge- invariant, the question is how one can possibly have a fermi 
vacuum specified by a fixed set of n F J s and still not violate gauge-invariance. The answer, as we 
shall see below, is that there is no fixed fermi vacuum in our theory with a given set of n F J s and 
the total wavefunction for the full theory (for example for the full ground state) is given by a 
sum J2{n a \ Wf) ® 4>n a ({A a }) where (fi n °- (A a ) is the (lowest energy) eigenfunction of the effective 
hamiltonian for gauge fields for given n F 's. 

We start by observing that |{n^}) is a wavefunction of the fermions. To construct wavefunctions 
of the full theory we need to consider tensor products of this with wavefunctions of the gauge fields. 
Denoting the latter by </>{n a }({A a }) we may write for the full wavefunction 

* nF (V;A) = IK-» ®0{«s.}({A«}) (HI) 

Since \n F ) is a simultaneous eigenstate of the fermionic operators OJJ = Tr(7 5 ^V ; n a ) an d since the 
O^'s commute with the full hamiltonian H, ^ the above wavefunction would also be an eigenstate 
of H (albeit a non-gauge-invariant one), provided <j){ n a\ satisfies the equation 



^-2§I + ^ ({U) 



<M{A a }) = E<t> n% ({Aa}), ia = K + n%\ Q (112) 



where H r F 9 is obtained from Hp in (|104| ) by substituting V reg {\) and K reg (X, A') repsectively for 
V(X) and K(X, A') in it. 

4 These operators commute with the quadratic part of Hf for massless quarks. It can be shown that their 
commutator with the quartic fermion term is down by a factor of 1/iV. 
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In the above equation the values of A a are in the iV-torus (A a £ [0, Ao))- Since this is a compact 
region we must supplement the differential equation (|112J ) by boundary conditions specifying what 
happens when any one of the eigenvalues (say A;,) reaches the value Ao- We will now argue that the 
appropriate boundary conditions are that both 4>{ n °-y and its gradient <9\ c 0{n a } are continuous: 

4>{n a F }{^b = A ) = (f){ n a F+S ab}(\ b = 0) 

(llo) 

9\ c <P{n a F }(h = Ao) = d Xc 4> {n a F+S a b} (X b = 0), c= 1,2, ...,b,...,N 

Equation ( jH3| ) is basically a consequence of invariance of the theory under topologically non-trivial 
gauge transformations and is connected with the phenomenon of 'spectral flow' ||. The point is 



best illustrated by considering the spectrum, eqn. (85), of the dirac operator for each colour a, 
which we reproduce here for convenience 

^ Q = ^sgn(a)(n+^) (114) 

Consider a one-parameter deformation 

H^X b (n)=fiXo, /i£[0, 1] (115) 

such that 

A 6 (0) = 0, A 6 (l) = A (116) 

Under this flow a right-handed (left-handed) energy level continously increases (decreases) as /i 
increases from towards 1. As ii reaches 1, however, the spectrum becomes identical to that at 
ii = (namely, (2-7r/L)x integers). This is in accord with the fact that \ b = and A& = Ao represent 
the same physical point on the A^-circle. We may wish to ask if the physics is always invariant 
under the change (At, = ())—► (\ b = Ao). In the pure Yang- Mills theory this is in fact what happens 
since A& = Ao can be gauge-transformed to X b = 0. In the presence of dynamical fermions, however, 
such a gauge transformation also changes the fermi momenta (see (p7|)) and in particular the fermi 
levels ((|90D,(pl|)). In this case the statement of invariance translates to the fact that physics should 
be invariant with respect to the change (Ab = 0, n b F ) —>■ (X b = Ao, n b F — 1). Indeed every term in Hp 
remains invariant under this change because in each of them A a and n F appear in the combination 



£a = A a + n*pAo for all o. The boundary conditions ( 113 ) reflect the fact that the amplitudes for 
the gauge field configuration (X b = Ao) at some given value of n b F is the same as that of (A& = 0) 
at n b F + 1. 

Note that since the differential equation (|112j ) is actually in terms of £ a = A a + n^Ao we can 
write 

<t>{n F }({K}) = u({X a + n a F }) (117) 

where u({^ a }) is a function of N real variables £ a £ (— oo, oo), satisfying the differential equation 



£-^+^ e9 ({a) 



The boundary conditions ( |113j ) basically ensure the continuity and smoothness of the function 
u ({£a}) a t integer values of £ a (it is enough to ensure continuity of the function and of the first 
derivative since the differential equation ( |118| ) ensures the continuity of the higher derivatives). It 
is indeed the boundary conditions ( |113|) that encode the information that (27rA a /Ao)'s are 'angles' 
and there is nothing special about the points (2-7rA a /Ao) = (or = 2tt) to warrant any discontinuity. 
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Let us now go back to the ansatz for the full wavefunction ( |111| ) of the theory. The gauge trans- 
formation of the wavefunction under f2({m a }) can be read off by making the inverse transformation 
W on the 'coordinates' tp a and A a . Using the fact that (t/j\Q\{rip}} = {ip\{n F — m a }}, we have 

vEv(<M) -» *tOM) = * nF (^ nt ,Af) = \{n a F-rn a })®4> {n a } {{\ a -m a }) (119) 



Using the boundary conditions ( |113| ) we can rewrite ( |119j ) as 

9 - * n (V,A) = ttty nt ,A nt ) = |{n| -m Q }} ® {n a_ ma} ({A o }) (120) 



The last equation not only tells us that the state (111) is not gauge- invariant but it also immediately 



tells us what the gauge invariant state is. It is simply 

l*)= E lu4}>^{n«}({A4) (121) 

To see the implication of the above, let us calculate the expectation value of the full hamiltonian 
H in the state \9). 

(^|^|^)=E { „^ } na/o A °^a^ n . } ({A a })F cff {n . } ({A a }) 

= II« IZo dtaU*({Z a }) (Za "^7^ + flP*(tf fl }) tl({U) 



where H e g is given by the expression in square brackets on the left hand side of ( [LI 2D . The last line 
shows that the presence of fermions forces a decompactification of the eigenvalues A a to the gauge- 
invariant combination A a + n F = £ a which is a real number £ (— oo, oo). It also tells us that there 
is no fixed fermi level in a compact gauge theory with dynamical gauge fields; rather, the effective 
theory of the gauge fields {A a } in the fermi vacuum is given by a density matrix p(n F , \ a \n'p, A^). 
The trace wih respect to the density matrices in the above calculations reduces to a single sum over 
{n F } because quantities like total energy etc. in the large A?" limit are diagonal in the occupation 
number representation. 

Connection with axial anomaly: 

The above discussion has an intimate connection with the existence of the well-known nonabelian 

axial anomaly in the present theory |13| 

£>%,5 = ^F^ (123) 

where J" 5 = tp 11 ^^^^ is the axial colour current. Let us denote by Q% the diagonal axial charges 
Jo ^0 5 = Jo '4 ) ^ a 1^4 ,a - Integrating both sides of ( |123j ) over all space and restricting to diagonal 
colour components we get (in the gauge d x A\ = 0, Af = \ a 5 a b) 

L 



dtQl + ig ! dx [A , Jo,5] aa = - I dx Fft (124) 

Jo vr Jo 



The commutator term on the left hand side vanishes in the fermi vacuum, while the integral on 
the right hand side is just the constant colour-diagonal component of the transverse electric field. 



With these inputs we may rewrite ( 124 ) as 



IE 
dtQl = -T^ ( 125 ) 
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where we have denoted (1/L) J dxFffi by E a . 

Let us now see how ( |125| ) arises from our discussion of the previous section. In the fermi vacuum, 
we have 

Q l\ Qa % Q% a _ oo ^ 

Qr — 2Jn=-oo 1) Ql — 2Jn=n£+l 1 

Using the exponential regularization once again, we get 

QW) = ~ + \ + n a F + ^ + o(6) = - + \ + f + o(e) (127) 

e z aq e z aq 

and 

Ql(e) = I-±- „■,-£ + o(e) =l-l-k + o(e ) (128) 

giving 

Ql = ^ + 1 (129) 

Thus we see that the regularized diagonal axial charges, QQ, depend on the gauge fields A a . In 
fact, as A a is tuned from to Ao (eqs. (|115| ) and ( |116| ) ), Q% changes by 2. The reason for 
this is the spectral flow which causes an existing left-handed energy level to be shifted to the 
right-handed energy levels. As a consequence of this spectral flow, the regularized diagonal axial 
charges are neither integers nor conserved. The best way to see the presence of anomaly (that is, 
nonconservation of Qg) is to work in terms of the extended variable £ a (see the example in ( |122j )) 
and to work out the Heisenberg operator equation of motion for the charges Q^ in the effective 
theory of the gauge fields. We get 



dtQl = -i{g/*N)d ia = 2E a /\ Q (131) 



Using ( J129D and ( pOD we get 

which is identical to ( |125|) .pl 

6. The effective theory of gauge degrees of freedom — description in terms of 
fermion bilocal operator 

We have seen that the effective theory of gauge degrees of freedom is described by a system of 
N nonrelativistic fermionsF] interacting with each other. The effective hamiltonian is given by 

H cS = -E a 2Tm d l + T Ea(WA + 1/2) 2 

~ 2LA*JV 2 ^a,6 [>P(Wab) + 1p(-W a b) + W ab {^' ' (w ab ) - tp' (-W a b)}\ , W ab = (^ a - ^)/A 

(132) 

5 In the last equality of (|13l[) we have identified E a with the operator (—i/N)d t . This is easy to see because on 
the wavefunctions (eq. (|ll7[)) the operator d\ a reduces to d^ a in the sense of eq. (122) 



6 The mapping of the l_j3laquette lattice gauge theory into free non-relativistic fermions is an old result JL4J . The 
method of proof used in |l4J can be trivially extended to any polygon approximation of the circle. In the continuum 
limit, as expected, the hamiltonian is given by the laplacian on the unitary group. The wave-functions are U(N) 
singlets in the absence of external sources. A recent discussion of YM theory on circle is given in [H]. 
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Introducing a new fermi field, ??(£,£), which satisfies the anticommutation relation 

m,t),vHC,t)} = 6(^-0, (133) 

we can build a field theoretic description for the effective theory: 

/OO ^OO /"OO 

d£ vHt,t) htffot) - d£ / dC'vHOviOK^vHeHO (134) 

-OO J —OO J —OO 



where 



;>2 , 5^0 / £ . Is2 



and 



'^-2l^^ + ¥<t + 5) (135> 

Av = K«-e'), K({) = ^_ J (i + esj)«<{) + *(-{» d36) 

The physical states of this system are required to satisfy the operator constraint 

coo 

d£ J&t) rifot) = N (137) 

which imposes the requirement of a fixed number (JV) of fermions. 



Because of the constraint ( 137 ) the physical states of this system are bosonic. A procedure 
for complete bosonization of systems like the present one has been developed in [To] , |16|| . The key 
ingredient in terms of which the bosonization is achieved is the bilocal operator $(£), defined by 

*«'(*) = V&t)vH?,t). (138) 

The bilocal operator satisfies a W^ algebra, which can be deduced from the fermion anticommu- 
tation relation (133): 



\b(t), *z&(tj\ = <5(6 - ei) %^(t) - <s& - £ 2 ) %&(*) (139) 

It also satisfies a quadratic constraint, which follows directly from ( |137| ): 

/OO 
df" %»(*) $?"H f ) = (iV + 1) %'(*), i.e. $ 2 = (N + 1) 6. (140) 

-OO 



Also in terms of the bilocal operator $ the constraint (137) reads 



Tr{\ - $) = N (141) 

where the symbol i Tr'' stands for the usual matrix trace, i.e. TrA = J^ d£Ae£. 

A description of the classical phase space of this system is obtained in terms of the expectation 
values ($) of the bilocal operator in Woo-coherent states. We shall denote these by <£. One can 
show that the classical phase space satisfies the constraints 

$ 2 = $ ) Tr(l - $) = N. (142) 



These constraints are classical analogues of the operator constraints ( [140D and (141) and define a 
coadjoint orbit of the Woo algebra (|139j ). One can now immediately write down a bosonic field 
theory action, following Kirillov's method of coadjoint orbits: 

if r r /"oo roo 

ScS= Nj dsdtTr ^( s ^)i d M s ^)^ d s^(s,t)})- J dtTr{h$) + J dt J dij d£'<$>t.i^tt' K tt> 

(143) 
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where 

h& = htS(£ - O (144) 

The integration region £ of (s, t) is the lower half plane, se(— oo, 0),ie(— oo, +00), and $>(s, t) satisfies 
the boundary conditions <&(s = 0, t) = 3>(£), &(s,t) = time-independent constant matrix. 

s— > — oo 

We note that the term quadratic in <J> in the action ( |143| ) breaks the residual part of W^ invariance 
that survives after introduction of the 'magnetic field' h. 

Because of the constraints (|142[), the allowed variations in $ are the W^ rotations: 



$ _> e iW§ e -iW &<$> + i[w, $] + ... (145) 

Making such a variation in 5 e g we deduce the classical equation of motion for 3>: 

/oo 
de^^^(K^,-K^) (146) 

-oo 



I 

N' 



One can easily check that this is identical to the operator equation obtained using ( |134j ) . Note that 
the factor of 1/N on the left hand side is there because ft = 1/N. 

The classical ground state of this system is described by such a time-independent solution of 
(|146| ) that satisfies the constraints (142). Excited states are described by arbitrary Woo rotations 



of the ground state. In general the description of these states is rather complicated. However, 



the semiclassical limit of a theory of the type described by the action (143) and constraints (142) 
admits of a simple fermi fluid description in the sense of a Hartree-Fock approximation. If we 
further restrict ourselves to small fluctuations of the fermi surface of a special type ('quadratic 
profile') then the fermi fluid description [|17]] further simplifies to a description in terms of the 
collective variables JLJI, the density p(£,t) and its conjugate momentum II(£, t). The semiclassical 
limit of the present theory is obtained when N — > oo. 



In the small L limit one can neglect the quartic fermion term in (134) and the theory considerably 
simplifies. We shall discuss this limit separately in a later section. 



7. The operator algebra of meson fluctuations 



In the previous sections we fixed our attention on the lowest energy state in the fermionic sector 
since fermionic fluctuations do not affect the effective theory of gauge fields in leading order in N. 
To obtain the meson spectrum in this theory, however, we need to consider fluctuations in the 
bilocal field. The allowed fluctuations are W"^ -rotations. Let us denote the expectation value of 
M{\) in the full vacuum of the theory by Mq(X) and the corresponding quantity for M(X) (see 
eqn. (fflj)) by Mq(X). Then, for fluctuations we must write 



M(X,t) = jW{\t)lVN M ( X)e -iW(X,t)/VN 

■ , , ( 147 ) 

= Mo(X) + ^[W(X, t),M (X)} + ^ [W(X, t), [W(X, t),M (A)]] + o(^=) 



Here A4(X) is defined as in (|7l|). By definition, A4q(X) satisfies the quadratic constraint .Mo (A) = 
.Mo (A). This implies the existence of the two projection operators P+ and P-, which are given by 

P + (X)=M (X), P_(A) = 1-M (A), P|(A) = P±(A) (148) 
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With the help of these projection operators we can distinguish four different kinds of operators that 
may be constructed out of A4(X,t): 

M±±(\,t) = P ± (X)M(X,t)P ± (X), M± T (X,t) = P ± (X)M(X,t)P T (X) (149) 

Using fli~47D and ([Uf) we get 

M ±T (X,t) = T -^W±T(X,t)+o(±) (150) 



M ++ (X,t) = Mo(X)-^W + „(X,t)W-+(X,t)+o(^=) (151) 

M--(X,t) = ±W-+(X,t)W + „(X,t) + o(j±=) (152) 

The W^c algebra satisfied by Ai( A), eqn. (|7q), implies an algebra for the operators M.±± (A), -M± T (A). 
It turns out that the first two of these generate a subalgebra of W^ , which acts on the other two. 
We give the full algebra below. 



M 



ia,jf3 
±±xj/ 



(A), Mit^'m 



fr* 



±xy 



JVl ±±xy 



(A), 



M 



JVl ±Tx'v l i A 



ia,j(3 



\'){pt'fWM^,(X) 

"Sp(X-X f )PSf(X)M^,(X) 



(A), M£±5 (A') 



A 

iV 

N 



y'WM^X(X)} 

(153) 

(154) 



Sp(X 



X')P^ \X)M^f y (X) 



(155) 



tioi,jf3 



m i :%(x) 



ii'cc'j'P 



M- + X ( A ') 



Ao 



N 



Sp(X 



^{p^fw^Zw 



-XV 



•-xy 



p*a>3 
Jr +xy' 



^(A)-M^(A)} 

(156) 

The commutators of all other combinations vanish. The first of the above is the statement that 
M±± (A) satisfy a subalgebra of W^. The second and third show that this algebra acts on M±zf (A). 
The last is best understood in terms of the fluctuations W± T (X). Let us substitute (|150|) — (|152j ) 
in the above algebra, retaining only the first non-trivial term in the fluctuations in —p= expansion. 

We see that (j\4. ++ {X) — M.q{X)) and M. (A), which are realized to this order in -4= in terms of 

the fluctuations W^ (A) and W ^(A), satisfy a subalgebra of W^ which acts on the fluctuations. 

The origin of this action is the commutation relation in (|156|) , which gives for the fluctuations 



Wf_^(A), W^'(X') = X 5 P (X-X')P^f(X)P 



\-xy' 



(A) 



(157) 



We see that a c- number term has appeared in the algebra of fluctuations ( |157| ) . Inn a correspond- 
ing phenomenon had occurred for meson fluctuations where the c-number term was identified with 
the central term of a Heisenberg algebra for the fluctuations. 



The c-number term in (157) depends on the vacuum expectation value of A4(X). We close 
this section by outlining briefly how this can be obtained. First we note that (by equation (fn])) 
A^o(A) = Mo(<7 _1 (A)), where g is the transformation that takes the eigenvalue distribution in 
the ground state of the effective gauge theory to the fiducial distribution (see eqn. (|65|)). Note 
that the decompactification of the eigenvalues discussed in section 5 does not affect this since 
the transformation g only sees the fractional parts of £ a 's (or A a 's). Now, g(X) can be expressed 
in terms of the density U(X) as in eqn. (|68|), It follows, therefore, that a computation of M.q(X) 
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involves knowing Mo (A) and Uq(X), where the latter is the vacuum expectation value in the effective 
gauge theory. Actually, both M (A) and Wo (A) can be obtained from a knowledge of the vacuum 
expectation value of a more general density variable p(£, t) defined on the entire real line: 

p(£,*) = -^ £<*(£- £«(<)) ( 158 ) 

a 

It is related to U(X, t) by 

oo 

U(X,t)= J2 p(A + mA ,t) (159) 

m=— oo 

Also, M xy (X, t), given by ( |99| ) and (100), can be rewritten in terms of p(£,t) as 



M (A, t ) = - T 5 — ^— piX,t 'ff V lr (160) 

where 

oo 

p(X,t;x-y)= J2 e^ m ^~y^ L p(X + mX ,t) (161) 



m=— oo 



Thus a knowledge of po(£) is sufficient to compute Mo (A) and Uo(X) and hence .Mo(A). A compu- 
tation of the value of /?o(£) i n the small L limit is given in the next section. 

8. The small L limit 



Here we are interested in the limit L — > 0. In this limit Ao = 2ir/gL — ► oo, i.e. as the length of 
the space circle becomes vanishingly small, the length of the internal circle becomes infinitely large. 
We then expect the dynamics to reside only in the gauge field degrees of freedom, since the fermionic 
fluctuations are further suppressed by factors of 1/L in addition to the 1/N suppression already 
present at large N. Now, the dynamics of the gauge fields is given by the effective hamiltonian 
H e g ( |132[) . The first two terms in H e g are of order 1/L. The last term is subleading in L and does 
not contribute. This is because its coefficient vanishes linearly as L — > while the quantity in the 
square brackets remains finite in this limitQ. Therefore to leading order in the L — > limit, the 
effective gauge field dynamics is described by a set of ./V free non-relativistic fermions moving on a 
line in a harmonic oscillator potential. This problem may be solved by the bosonization methods 
introduced and extensively discussed in |l(| [l6|] in the closely related context of the c = 1 matrix 
model, in which the only difference is that the potential has the wrong sign. 

In Sec. 6 we have introduced this formalism of bosonization for the more general problem 
of interacting fermions. In the small L limit the quartic fermion term in H., eqn. ( |134| ), can be 
neglected and the theory considerably simplifies. The relevant equation of motion in this limit is 

±d t $(t) = Mt),h}. (162) 

7 This can be seen by using the property of the function ip{x), ip(x) — ip(x + 1) — 1/x. At small x ip( x ) behaves as 
— 1/x since ip(l) is finite. Using this, it follows that the quantity in the square brakcets in H g vanishes as w a b — * 0. 
It is this last limit that is of relevance here since the maximum separation between two of the scaled eigenvalues 
£ a /Ao and £f,/Ao goes as \[L. 
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The classical ground state of this system, $0, is described by a time- independent solution of ( 162 ) 
which satisfies the constraints 

<$>l = $ 0> Tr{\ -Q )=N (163) 

Let us parametrize such a solution as 

oo 

$o,«' = 6(Z - - E CnM0<PUO (164) 

n=0 



This satisfies the equation of motion ( |162| ) provided the </> n 's satisfy the Schrodinger equation for 
a harmonic oscillator: 

n\ n fi 1 -1 

E n <j>n (165) 






2LiV 2 ? iV v A n 2 



y n — J - J n K rn 



Assuming that </> n 's are orthonormal and imposing the constraints (|163|) on (164), we get 



c 2 n = en, J2c n = N (166) 

n=0 

The first of ( |166j ) says that any given level is either occupied or unoccupied (fermions!), while 
the second says that precisely N levels are occupied. Since the energy of the level labelled by n 
goes as n, minimising energy in the ground state implies that the first N levels are occupied, i.e., 

c n = e(-n + N-l). (167) 

The density p(£, t) is given by the formula 

p(ti,t) = ±(l-m\c (168) 



H N-l 



For the ground state we see that 

M6 = ^ E MOfniO- (169) 

n=0 

Excited states (glueballs) are described by W m rotations of the ground state constructed above. 
Writing 

$(£) = e iA W$ oe - iA W (170) 

and neglecting quadratic and higher terms in A(i), we see that A(i) satisfies the equation of motion 

[±dtA(t)-[A(t),h],$ o ] = (171) 

This leads to 

^d t A ±T (t) = [A ±T (t),h] (172) 

where 

A + _(t) = <DoA(t)(l-$o), A_+(f) = (1 - *o)A(t)* (173) 

This construction is identical to the one just described for meson fluctuations and everything 
discussed there also applies here. Defining 

/OO /"OO 

dU ^«(0(^(%i(O (174) 

-oo J —oo 
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and using the ground state solution <J>o obtained above one can show that A"™(£) vanishes if n < N 
or m > N. Similarly, one can also show that A"™(£) vanishes if n > N or m < N. Moreover, these 
two satisfy the algebra 

[A™, A*?'] = 9(n - N)9{-m + N - l)<r ,m V m (175) 

Thus, for appropriate ranges of indices, A"™ and A!! 1 " form a canonically conjugate pair. 
The equation of motion for A^2 is 

^d t Al™(t) = -{E n - E m )A^{t) (176) 

which is trivially solved by 

Thus the spectrum of glueballs consists of 'discrete' states labelled by two integers (n,m) JT^]. A 
residual wedge subalgebra of the original W^ acts on this tower of states and moves them around. 
This subalgebra is generated by 

$++ = $ $$o, (178) 

and 

*__ = (l-$o)*(l-*o), (179) 

exactly like the situation discussed in the previous section. $"™ is nonvanishing only for n,m > N, 
while <3>™™ is so for n,m < N. Their action on A_| is given by 

[$ nm A rSm'} = q(j i _ ^n'm^nmf (180) 

[$^, A%™'\ = -6{-m' + N - l)5 nm ' A r l m (181) 

There is a similar action on A"™. Thus the spectrum in the limit L — ► consists of a tower of 
discrete states (glueballs) which realize a representation of a wedge subalgebra of W^ algebra. 

9. The large L limit 

In the limit L — v oo the pure Yang-Mills part of the hamiltonian behaves like a system of 
infinitely massive particles. Fluctuations in {A a } are therefore small, and so we can fix our attention 
on the ground state only in the Yang-Mills sector. Also as L — > oo, Ao = lix jgL — ► 0, i.e. the 
length of internal circle becomes vanishingly small. We then expect only the zero mode of M(A) 
to survive. 

Actually, in this limit M(A) is not a convenient variable to work with since it is not periodic in 
A. However, M(A), to which M(\) is related as in (J5J), is periodic and so this is the appropriate 
variable to work with in the large L limit. In terms of M(A) the equation of motion (78) reads 



Hdt-(l/L)J2aPad Xa ) M%M(\,t) 



j 5 id + m 7 °, M (A, t) + %- /;° d\' U(X',t) fn dz K(\' -\;x-z) (182) 

Ixy z L 

/ — — \ice,jB — / — — \ia,iP 

[M XZ (X', t)M zy (X, t)) - K(X - X'y - z) (M XZ (X, t)M zy {X' , t)) 
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where K(X;x) = e t9xX K(X;x),K(X;x) being given by (|80|). Moreover, M(X) satisfies constraints 
that are identical to (|57|) and (paj): 

M 2 (X) = M(X), (183) 

Tr(l - M(X)) = 0. (184) 

Let us now expand M(X, t) in modes in A: 

M(X,t)= J2 e~^~M^ n \t). (185) 



Retaining only its zero mode in the equation of motion (182) and in the constraints ( |183| ) and 
dH), we get 



j 5 id + rn^°,M^(t) +?- M(°\t),M'(°\t) 


(186) 


(m(°)) 2 = m(°), 


(187) 


Tr(l-M (0) ) = 0, 


(188) 



where M xy (t) = Mxy (t)\x — y\. These are precisely the equations of motion and the constraints 
that one gets in the gauge A\ = for the fermion bilocal for QCD2 on a plane. The nonzero modes 
of M(X) decouple because excitation of these modes makes the corresponding states infinitely heavy 
in the limit L — > 00. 

We mention here that in ||] we solved the model defined by the equations ( |186| ), ( |187j ) and 
(|188| ) in the manifestly Lorentz-invariant gauge A + = 0. We showed there that the spectrum of 
excitations of this model is identical to the 'tHooft spectrum for mesons and that these particles 
realize a representation of a wedge subalgebra of the Woo-algebra which the fermion bilocal fields 
satisfy. 

10. Concluding Remarks 



To conclude we remark that some of the ideas presented in this paper may well apply to finite 
temperature QCD. It is well-known that field theories at finite temperature (T = 1/P) can be 
mapped to Euclidean field theories with compactified time (with period (3). On the other hand, 
for Euclidean field theories, time is just any one of the dimensions. Therefore, if we consider the 
Euclidean version of our theory, then we can learn about the spectrum of glueballs and mesons at 
j3 — > and j3 — ► 00 from our analysis at L — > and L — > 00 respectively. Indeed this is perhaps 
the simplest model involving dynamical quarks and gauge fields where one can learn about their 
excitations and interactions at finite temperature in an analytic fashion. Implications of this model 
for QCD at finite temperature are being worked out. 

The second remark is that by tuning L in this theory we seem to achieve a 'melting' of topological 
degrees of freedom to propagating degrees of freedom. As we have seen, in the small L limit (/3 — ► 0) 
mesons are frozen, and the theory is described in terms of the quantum mechanical degrees of 
freedom (X a (t)) which do not depend on space. As we keep increasing L, the mesonic fluctuations 
come into play. It is a very interesting open question how the theory behaves at intermediate values 
of L. 
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The third remark is about what we can learn from this model vis-a-vis string theory. The 
interaction between glueballs and mesons in this theory can be compared with an interaction 
between closed and open strings respectively where the closed string sector depends only on time 
and not on space. The freezing of the propagating (open string) degrees of freedom at L — > 
tempts one to identify the remaining, purely quantum mechanical, closed string degrees of freedom 
as the degrees of freedom of a 'topological' string theory. There is also the interesting possibility 
that this theory has a three dimensional (x, t, £) description where £ is a new dimension arising out 
of the eigenvalues in the gauge sector. We leave these interesting issues for future work. 
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